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Abstract 

We calculate the spin and charge dynamical susceptibilities of a strongly correlated impurity 
model in a renormalised perturbation theory. The irreducible for vertices for the quasiparticle 
scattering are deduced from the renormalised parameters, which have been calculated by fitting 
of the low-lying levels of a numerical renormalization group (NRG) calculation to those of an 
effective Anderson model. The susceptibilities are asymptotically exact in the low frequency limit 
and satisfy the Korringa-Shiba relation. By comparing the results with those calculated from a 
direct NRG calculation, we show that the renormalised perturbation theory (RPT) description 
gives a very good description of spin dynamics for all values of the local interaction U, not only at 
low frequencies but over the whole relevant frequency range. 

In the presence of a magnetic field the approach can be generalised using field dependent renor- 
malized parameters to calculate the transverse and parallel spin dynamic susceptibilities. The RPT 
results give accurate results for the spin dynamics over the whole frequency range and for all values 
of the magnetic field. 
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I. INTRODUCTION 



The scattering of electrons with low energy spin excitations is responsible for the 
anomalous behaviour of many systems which are described as being strongly correlated. In 
magnetic impurity systems this scattering is enhanced leading to the Kondo effect and the 
many-body resonance (Kondo resonance) which develops at low temperatures at the Fermi 
level. It is also the origin of the large mass enhancements in heavy fermion systems, and 
the strong renormalisation of the energy bands at the Fermi level. At a quantum critical 
point, where the critical temperature for magnetic ordering is reduced to zero by pressure 
or alloying, it is the local low energy fluctuations which cause the breakdown of Fermi liquid 
theory. In the cuprate high temperature superconductors, it is the scattering with the spin 
fluctuations, after the antiferromagnetic order has been destroyed through doping, which is 
the likely cause of the anomalous electronic behaviour in the normal state. The exchange 
of the magnetic fluctuations between the conduction electrons may be the origin of the 
electron attraction leading to superconductivity. It is important, therefore, to understand 
the nature of these fluctuations in the absence of any long range magnetic order. Here, we 
study the spin dynamics in a strongly correlated local model, the impurity Anderson model 
Q], and calculate local spin dynamics using a renormalised perturbation theory expansion. 
We start by describing the model and the technical aspects of the approach in this section 
of the paper, and later present results for various parameter regimes, both with and without 
an applied magnetic field. 

The Anderson model has the form. 



where is the energy of the impurity level, U is the interaction at the impurity site, and Vk 
the hybridization matrix element to a band of conduction electrons with energy e^. In the 
wide band limit the hybridization weighted density of states, A{uj) = nJ2k — ^k), 

can be taken as a constant A. With this assumption the impurity one-electron with spin a 
Green's function is given by 



Ham = (^ddU^ + Und,^nd,i + 5](Vfc4cfc,a + K(i,ada) + ^k^acl^^Ck, 



(1) 





(2) 



where T,„{u) is the self-energy. 



We have demonstrated in earher work that the low energy fixed point of this model 
can be interpreted in terms of the same model but with renormalised parameters [J, ^ . The 
renormalised model takes the form, 

^AM = ^ddid^ + U : nd.trid.i : + I^(Vfe4cfc,a + Vk4,M + ^k,acl^Ck^„, (3) 

k,>7 k,cr 

where the colon brackets indicate that the expression within them must be normal-ordered. 
The three parameters which specify this renormalised model are e^, A and U, defined by 

Q = z(erf + S(0)), A = zA, U = z^T^^{0,0,0,0), (4) 

where z is given hj z = 1/(1 — S'(0)) and r^i{uJi, 002, (^3,004) is the local 4-vertex. jj]. The 
renormalised Hamiltonian Q describes only the low energy fixed point and the leading 
irrelevant corrections. A complete description of the original model can be obtained by 
working in the Lagrangian formulation. The Lagrangian of the original model CAui^d, A, U) 
can be rewritten in terms of the renormalised model plus counter-terms, 

£AM(ed, A, U) = CAui^d, A, U) + £ct(Ai, A2, A3), (5) 

where £ct(Ai, A2, A3) describes counter-terms. In this formulation a renormalised pertur- 
bation expansion can be carried out in terms of the quasiparticle propogaters in powers 
of the renormalised interaction U. The role of the three parameters Ai, A2 and A3 in the 



counter-term is to cancel any further renormalisation o: 



which are taken to be fully renormalised from the start 



m 



the parameters , e^, A and U, 



In terms of the renormalised parameters the exact results for the impurity spin and charge 
susceptibilities take a simple form, 

xs = ^p(o)(i + f/p(o)), xc = lp{m-um), (6) 

where Xs(0) is in units of {gji-QY, Xc(0) differs by a factor from the usual definition, so that 
for the non-interacting system, Xc(0) = Xs(0), and p{uj) is the local density of states for the 
non-interacting quasiparticles. 



These results correspond to a first order calculation in the renormalized perturbation 
expansion in powers of U 



We can choose to work solely in terms of the parameters, q, A and U, and determine 
them by comparison with experiments, as one does in a phenomenological Fermi liquid the- 
ory. To determine them in terms of the bare parameters of the orginal model using equation 
(0} it would seem that we need almost the complete solution of the model, including the 
dynamics. However, they can be calculated without recourse to these formulae by fitting 
the low energy excitations of the numerical renormalisation group calculations to those of 
the effective low energy model given in equation (jH)), as was done essentially in the original 
NRG calculations The results for q, A and U, for various parameter regimes of the 

model have been described in earlier work 



We wish to generalise the results for the spin and charge susceptibilities given in 
equation ^ to obtain the corresponding dynamic susceptibilities: the spin susceptibil- 
ity, Xs{^), which is given by the Fourier transform of the double-time Green's functions of 
either {{d\{t)d^{t); d\{t')d^{t'))) /2 (t ransverse susceptibility), or ((nd|(t) — ^^^(t); Tidt(^') " 
nd;(t')))/4 (paralllel susceptibility), as they are equal in the absence of a magnetic field, and 
the charge susceptibility, Xc{^^) = ((^dt(^) + "-di(^) : n^iit') + ^dj.(^')))/4- To calculate these 
we need to look at the processes involving repeated quasiparticle scattering. A typical dia- 
gram in the renormalised perturbation expansion corresponding to the repeated scattering 
of a spin up quasiparticle with a spin down quasihole is shown in figure Q These scattering 
terms at u; = have implicitly been taken into account in the results for the static sus- 
ceptibilities in ®, and so we must introduce the counter-term A3 to prevent overcounting. 



We define an irreducible interaction vertex U-, 



H 



U — X-i, and sum the series to derive an 



expression for the dynamic spin susceptibility, 



21 



(8) 



where O.SIIJ^Uu;) is the spin dynamic charge susceptibility for the non-interacting quasi- 



particles. Expressions for n^j(c(j) for the particle- hole symmetric model are given in 
the Appendix. The approach is similar to the random phase approximation (RPA) in 
terms of quasiparticles 0,0] and goes over to the RPA expression in the weak couphng limit. 




FIG. 1: The repeated scattering of a quasiparticle t and a quasihole J, which is taken into account 
in the RPT calculation of the dynamic spin susceptibility Xs{^)- 



At zero frequency we can identify the result with the spin susceptibility given in equation 

X.(0) = - ^j? = ^(1 + Um), (9) 
The effective particle-hole vertex can then be deduced in terms of U, 

''^ l + f/p(0) ^ ' 

In terms of conventional Landau Fermi liquid notation, p{0)IJf^^^ = —Fq. The resulting 
expression for Xs{^^) in terms of U is 

1 nl\{u){i + um) .... 



2 1 + Up{0) - UIl 



This result is asymptotically exact as — ^ as Xs{^) p(0)(l + f^p(0))/2 (by construction), 
and its imaginary part satisfies the Korringa-Shiba relation, 

limi^^^ = ^(l + £/«0)f = 2x2(0), (12) 



which is an exact relation for the Anderson model lllll. We note that the result (IllD . in the 

n 

large U limit where Up{0) ^ 1, is equivalent to that derived by Kuramoto and Miyake |12| ) 
based on a semi-phenomenological Fermi liquid theory. 

We repeat this exercise but this time consider the repeated scattering two opposite spin 
quasiparticles, as illustrated in figure |21 and define a corresponding irreducible interaction 
vertex ?7p|. The sum of this series gives for the charge dynamic susceptibility. 



= — T-^r — ■ (13) 




FIG. 2: The repeated scattering of a quasiparticle t and a quasiparticle | which is summed in the 
RPT calculation of the dynamic charge susceptibility Xc{^^)- 

Identifying the u = the result with the charge susceptibihty in equation (jH)), 

XM = i = M(i _ u~p^o)). (14) 

The effective particle-particle vertex can then be deduced in terms of U, 



U_ 

1 -"^-(0) 



Kl = -r^^ (15) 



In terms of conventional Landau Fermi liquid notation, p{0)U^J_= Fq, and the values of 
and Fq are such as to satisfy the forward scattering sum rule [10|, 

° +7^ = 0- (16) 



1 + F^ 1 + F^ 

In figure El we give the two irreducible vertices, f/^j/yrA and t/p|/7rA, for the symmetric 
Anderson model as a function of fZ/vrA. It can be seen that there is a rapid divergence 
between the interactions in the spin and charge channels with increasing U. The vertex 
fj^l in the spin channel decreases monotonically with increase of U and in the limit of 
large U it tends to the value 2Tk, where Tk is the Kondo temperature given by Tk = 



([/A/2)e~'^^/'^^"'"'^^/^^. The vertex JJ^j in the charge channel, on the other hand, increases 
monotonically though remains finite and diverges only in the limit U oo. 
The RPT result for the dynamic charge susceptibility is 

1 iiaM(i-[/p(o)) 

Yc = %z ^^^-f , 17 

^ ' 2 i-f/p(o) + f/nf^H' ^ ' 

where np|(u;) for the symmetric model is given in the Appendix. This result is again asymp- 
totically exact and satisfies the relevant Korringa-Shiba relation. In the next sections we 
evaluate these expressions for the particle-hole symmetric model, and compare the results 
with those calculated directly from numerical renormalisation group calculations. 




FIG. 3: The irreducible vertices, Uf^j^{h) (full line) and U^j^ (h) (dashed line) as a function of U /vrA 
for the symmetric Anderson model. 

II. DYNAMIC SUSCEPTIBILITIES IN THE ABSENCE OF A MAGNETIC 
FIELD 

We have shown that the RPT expressions for Xsi^^) and Xc{^) a-i's asymptotically exact 
in the limit — 0, but we would like to know over what range of finite u they give accurate 
results. To estimate this we compare them with the corresponding results derived from a 
direct numerical renormalisation group (NRG) calculation for the dynamics [l^l. 11^. These 
spectra can be calculated using the same methods which have been developed for calculating 
the spectral density of the one-electron impurity Green's function. These direct calculations 
give a good overall picture of the behaviour over all energy scales, but have some inaccuracies 
stemming mainly from the use of a discrete spectrum for the conduction electrons, and the 
logarithmic broadening which is used to give the final continuous spectrum. To gauge the 
accuracy of the direct NRG calculations we calculate Xs{^) using this technique for the 
non-interacting system {U = 0) and which we can then compare with the exact result which 
for the particle-hole symmetric model corresponds to equations and p7|) with U = 0, 
Q = and A = 0. We take A = 0.015708 for these and all subsequent calculations. 

The results for imaginary part of the retarded susceptibility, ImXsi^^) are shown in figure 
m We see that the direct NRG results give a reasonably good representation of the overall 
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FIG. 4: A comparison of the results of an NRG calculation (dashed line) of the imaginary part of 
the spin susceptibility Xsl"^) for the non-interacting symmetric Anderson model (U = 0) compared 
with an evaluation of the exact result (full line). 

behaviour. However, there is some overestimate of the gradient at = 0, leading to a 10% 
error in the Korringa-Shiba relation, and about a 5% overestimate of the peak height. 

We can now make a comparison of the direct NRG results in the interacting case with 
the corresponding RPT results. We look at a strong correlation situation first of all with 
U/tiA = 3.0, in which the low energy charge excitations are almost competely suppressed. 
The values of If and A are calculated from the NRG levels as described in earlier work. 
The results for lmxs{^) are shown in figure |S1 The peak structure is very much enhanced 
compared with the non-interacting case and occurs on a scale ~ A = AT-k/tc, where Tk is 
the Kondo temperature. There is very good agreement between the two sets of results over 
the whole range of this low-lying peak structure. It would appear that the RPT results may 
be in fact more accurate than the direct NRG results, as they satisfy the Korringa-Shiba 
relation exactly and, as we saw from the results in figure |31 the direct NRG results tend 
to overestimate the peak heights. The results for the real part of Xs{^) is shown in figure 
inland there is good agreement between the results but the direct NRG underestimates the 
value of Xs(0) by about 4%. The corresponding dynamic charge susceptibilities are very 
small over this range. 




FIG. 5: Results for the imaginary part of the dynamic spin susceptibihty Xs{'^)/'^ for the symmetric 
model in the Kondo regime U /vrA = 3.0. The dashed curve corresponds to the NRG results and 
the full line to the results using the RPT expression, equation 



We take a smaller value of U such that Uj-nl^ = 1.0 to see spin and charge peaks of 
the same order. The NRG and RPT results for both ImXs(c(j) and Imxd^) are shown in 
figure [7| The enhancement of the spin peaks are down by a factor of 10 from the strong 
correlation case shown in figure and the u range shown is greater by the same factor. 
In this moderately correlated case the peaks in the charge spectrum are reduced but can 
be seen on the same scale as the spin peaks. Again there is very good agreement between 
the NRG and RPT results over this whole range, and the evidence from the earlier results 
suggest that the RPT results may be the more acccurate ones over this range. 



We take a case of weak correlation with U /vrA = 0.5. In this case it should be meaningful 
to make a comparision with the standard RPA results, which correspond to taking the bare 
parameters in the RPT formulae; = U, A = A. The results of all three calculations 
for the imaginary part of the dynamic spin susceptibility are shown in figure |H1 AH three 




FIG. 6: Results for the real part of the dynamic spin susceptibility Xsii^) for the symmetric model 
in the Kondo regime U /vrA = 3.0. The dashed curve corresponds to the NRG results and the full 
line to the results using the RPT expression, equation 




FIG. 7: Results for the imaginary parts of both the dynamic spin and charge susceptibilities sus- 
ceptibility Xs{^) foi" the symmetric model with U /vrA = 1.0. The dashed curve (spin) and the dot- 
dashed curve (charge) corresponds to the NRG results and the full lines (thick/spin:thin/charge) 
to the results using the RPT expressions. 
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FIG. 8: The imaginary part of the dynamic spin susceptibihty for a weak correlation case of the 
symmetic model with [//vrA = 0.5. The three sets of results are from (i) NRG calculations (dashed 
line), (ii) RPT (full line) and (iii) RPA (dot-dashed hue). 

curves are qualitatively similar, and there is very good agreement between them regarding 
the peak positions. The differences between the NRG and RPT are similar to those noted 
in figures 4 and 6. If we compare the peaks with those in the non-interacting case shown in 
figure 01 we see that they are enhanced only by a factor of the order 1.5. The peak heights 
are overestimated in the RPA results. 

In figure M we present the results for ImXs(ct;)/7rx^(0), for the four sets of results 
U /ttA = 0.0, 0.5, 1.0, 3.0. All sets have the value 2 aX u = corresponding to the Korringa- 
Shiba relation. With increasing U there is the same narrowing to a scale set by the Kondo 
temperature as in the one-electron spectral density of the impurity Green's function. Com- 
bining the the Kramers-Kronig relation, J^^lmx{uj) /uj du = 7rRex(0), with the Korringa- 
Shiba relation gives 



1 



oo 



Im XsA^) 



duj 



1 27rA 



(18) 



— oo 



UJ 



Xs,M (l±t>/7rA)' 



The total weight under the peak in figure IHl therefore, when U/ttA ^ 1 and U /ttA — > 1, is 
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FIG. 9: Results for ImXs('^)/7i"Xs(0) calculated using the RPT for the symmetric model and a 
range of values of U. 

equal to 4Tk. The sum rule (II 8j) is satisfied precisely in the results shown. 

These results contrast to the corresponding results for the dynamic charge susceptibility 
shown in figure ^1 Again due to the Korringa-Shiba relation all sets have the value of 
2 at u; = 0. The suppression of the charge fluctuations in this case is indicated by the 
development of symmetric high energy peaks, with values very much greater than the value 
at u; = 0. From equation (fTHj) it can be seen that as f/ — > oo and If /ttA 1 that the 
area under the curves diverges in this limit. However, though qualitatively correct the RPT 
restricted to this class of diagrams cannot accurately describe the behaviour on this high 
energy scale. 



III. DYNAMIC SUSCEPTIBILITIES IN A MAGNETIC FIELD 

We now consider the dynamic spin response for the symmetric model in the presence of an 
applied magnetic field H. The low energy behaviour can still be described by a renormalised 
Anderson model with parameters which now depend on the magnetic field |15i]. The field 
dependent parameters h, A{h) and U{h), are defined by 

h = z{h){-^ + j:{0,h)), A{h)=z{h)A, U{h) = z^{h)T^^{0, 0,0,0 : h), (19) 

where z{h) is given by z{h) = 1/(1 — S'(0,/i)) and h = qunH/2. These parameters can 
again be determined from the levels of a NRG calculation |l5L llq . The induced impurity 
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FIG. 10: Results for Imxc(w)/7rx^(0) calculated using the RPT for the symmetric model and the 
same range of values of U as for the spin susceptibility results shown in figure IHl 



magnetisation Md{h) = {g^BY^{h) at T = 0, where m{h) is given exactly by 

mih) = -tan"^ | | , 



(20) 



a result which follows from the Friedel sum rule \V^. The renormalized parameters, fj{h) = 
h/h, A{h) and U{h), in the presence of a magnetic field H, relative to their bare values, for 
U/nA = 3.0 are given in figure ITTl The parameters are not al independent and U{h) can 
be deduced from the other two parameters from the equation 3| , 

dfj{h) hfi{h) dA{h) 



1 + U{h)p{0,h) = fi{h)+h- 



(21) 



dh A{h) dh ' 

where p{uj, h) is the quasiparticle density of states in the presence of the magnetic field and 
is given by 

P>.ft)- , (22) 



(cU-/l(/l))2 + A2(/l) 

As the magnetic field increases the parameters change from the Kondo regime {h < Tk), 
where fjQi) = 2, and U{h)p{0, h) = 1, through a local moment regime {h > Tk), where fj{h) 
increases with h, and U {h) is enhanced to values greater than the bare value U, and finally 
to the almost free model when h > U, and the parameters revert to their bare values. 



With a magnetic field present we have to distinguish between the parallel and transverse 
spin susceptibilities, Xs,\\{^^, h) and Xs,i.{^^ h). The static values at T = are given exactly 



1 



FIG. 11: The field dependence of the ratio of the renormalised parameters relative to their bare 
values, r]{h) = h/h, A{h)/A and U{h)/U for CZ/vrA = 3.0 plotted as a function of ln(/i/TK). 



by 



Xs,||(0, h) = ^p(0, h){l + U{h)p{0, h)), x.,±(0, h) = '^ 



1 



-tan- m] 

27r/i \A{h) J 



(23) 



To calculate the frequency dependence of these spin response functions from the quasipar- 
ticle scattering we need to calculate the correponding irreducible particle-hole vertices again 
as they will now depend on the value of the magnetic field. The transverse susceptibility 
corresponding to a excitation with a change in the z-component of the impurity spin of 1 
involves the scattering of a quasiparticle with spin up with a quasihole with spin down, and 
the irreducible vertex in this channel we denote by t/^|(/i). We derive its value using the 
same type of argument as earlier by requiring the the sum of the series for uj = 0, should 
equal the transverse static susceptibility given in equation ((221) • We obtain the result, 

nh{fi{h) - 1) 



tan-\h{h)/A{h))' 



(24) 



Using the results, 
lim{fj{h) 



1) = Um = lim 



tan-\h{h)/A{h)) 
h 



1 + um, 



(25) 



it can be shown that Uj^Uh), goes over to the earlier result in equation (fTUI) in the limit h ^ 0. 



The resulting expression for Xs,L{i^,h) is as given in equation ()11|) with the /i-dependence 
included in both the scattering vertex ijf}{h) and H' 



ln(h/r^) 



FIG. 12: The irreducible vertices, Uj^]^{h) (full line), Uj^J{h) (dot-dashed line), and UpJ^{h) (dashed 
line) as a function of In (h/Tx) for U/ttA = 3.0. 

The parallel dynamic susceptibility corresponds to repeated scattering of a quasiparticle 
with spin a and a quasihole with the same spin component a, and hence no change in 
the z-component of the spin. We denote the corresponding irreducible vertex by Uj^j{h), 
corresponding to cr =|, but for the symmetric model it is the same as for a =|. As the 
static parallel susceptibility has the same form as earlier in equation ©, apart for the h 
dependence, the result is 

Ufr(.l^ih) = , (26) 

with the corresponding result for the particle-particle vertex. Using the results from figure 
fTT| we can deduce the three types of vertices as a function of /i/vrA for U/tx/S. = 3.0. The 
results are shown in figure d Initially for small h, Ul\{h) = Ul]{h) = 2Tk, and then slowly 
diverge with Ulj{h) > Uj^j{h), and then slowly diverge with but converge again for h > U to 
the bare value but converge again ioi h> U to the bare value U . The charge susceptibility, 
which is initially very large, decreases and converges to the bare value U for h > U . 

The transverse dynamic spin susceptibility is antisymmetric in the absence of a field, 
ImXs^_L(ct;, 0) = — Imxs,±(— cu, 0), as can be seen in the results in figures EHHl and as a conse- 
quence the total integrated spectral weight is zero. This is no longer the case when h 0, 
and the total spectral weight w±{h) is equal to the average z-component of the impurity spin 
m{h) = 0.5{nii-\ — ridi). At T = the RPT result for h) satisfies this result exactly 
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FIG. 13: The NRG results for the imaginary parts of Xs,Li'-^,h) (dashed hne) and (dot- 
dashed) as a function of uj for U = compared with the exact results (_L-thick full line, ||-thin full 
line) for h = 0.47rA. 



as 

w±{h) = lim uJXs,±{^,h) = 0.5 lim ujI[^\{uj , h) = — tan^^ (a) ' ^'^'^^ 

which is equal to m{h) from equation ()20|). 

The total spectral weight w\\{h) for the parallel susceptibility Xs,\\{^, h) on the other hand 
vanishes as 

W||(/i) = \hn^uJXs,\\iuJ, h) = 0.5 jim^a;n^|(cj, h) = 0. (28) 

To gauge the errors in NRG results [l^ due to the use of a discrete spectrum and log- 
arithmic broadening we again compare them with exact results for U = 0. The results of 
this comparison for h = 0.4A are shown in figure ^1 We see that the NRG results give 
the correct overall forms for the imaginary parts for both the transverse and parallel sus- 
ceptibilities. However, due to the discrete spectrum and broadening on this relatively high 
energy scale, there is there is some deviation from the exact results. The asymmetry of the 
results for the transverse case due to the field is apparent. In figure IT^ we compare the NRG 
results for the imaginary part of the transverse susceptibility with the exact results for an 
extremely large magnetic field h = 5 A. On this scale the impurity is effectively an almost 
localised and nearly fully polarised local moment. The spectrum, therefore, consists of a 
single peak at the energy required to flip the spin in the applied field, Up ~ 2h. Whereas 
the exact result has a sharp peak the NRG is much broader on this high energy scale, but 
the peak position and total spectral weight of the NRG result are correct. 
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FIG. 14: The NRG results for the imaginary parts of Xs,Li^, h) (dashed hne) and as a function of 
uj for U = compared with the exact results (full line) for h = SvrA. 
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FIG. 15: The RPT results for the imaginary parts of Xs,L{^^,h) (thick full line) and Xs,||(^^7^) 
(thin full line) for h = 0.4Tk as a function of uj for U /vrA = 3.0 compared with the NRG results 
(_L-dashed line, ||-dot-dashed line). 

We now compare the NRG results with those of the RPT results in the strong coupling 
case corresponding to the Kondo limit. In figure IT5l we compare the imaginary parts of the 
parallel and transverse spin susceptibilities for a magnetic field h = 0.4Tk. We see that 
there is good agreement between the two sets of results. Comparing the differences with 
those seen in figure IT^ where the NRG results are compared with the exact ones, one would 
conclude that the RPT results are probably the more accurate ones over this range. 

The question arises as to whether the RPT results can give similarly good results for 
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FIG. 16: The RPT results for the imaginary parts of Xs,±{^Th) (thick full line) for h = 2Tk as a 
function of to for U/nA = 3.0 compared with the NRG results (dashed line). 

higher magnetic field values. As we increase the magnetic field the Kondo effect is suppressed 
and for extremely large fields the renormalised parameters revert to their bare values as can 
be seen in the results shown in figure ^2 We increase the applied field to test the results 
over this range of fields. We concentrate on the transverse susceptibility as this can change 
significantly with magnetic field as we pass from the Kondo regime to the local moment 
regime, while the parallel susceptibility mainly gets suppressed with increase of field. 

In figurelTHlwe show the imaginary part of the transverse susceptibility for U / ttA = 3.0 for 
h = 2Tk (Iu/i/Tk = 0.693). There is quite a dramatic change from the results for h = 0.4Tk, 
shown in figure E| The spectral weight in the region < is almost completely suppressed, 
and there is a single sharp peak in the spectrum for u > 0. Again there is remarkable 
agreement of the RPT results with those obtained from the direct NRG calculations. 

For values of the magnetic field such that ln{h/TK) ~ 3.0, we see from figure ITTl that 
U{h) is enhanced over its bare value U, corresponding to the local moment regime. It is 
interesting to see whether the Fermi liquid theory can be applied also in this regime. In 
figure El we plot the transverse spin spectral density calculated with the NRG and RPT 
for /i/Tk = 20 (Iu/i/Tk = 3.0). We see that remarkably that there is good agreement in 
this regime as well, both the peak positions and total spectral weight are in agreement. 
The peak position corresponds to 0.86 x 2h Though it can be seen from figure ITTl that the 
renormalised field parameter h{h) is enhanced over h, the peak position is not at 2h, and is 
somewhat less than the simple Zeeman splitting 2h. The NRG results are much broader but 
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FIG. 17: The RPT results for the imaginary parts of Xs,A.{^, h) (thick full line) for h = 20Tk as a 
function of uj for U/t\ 
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FIG. 18: The RPT results for the imaginary parts of Xs,±('^)^) (thick full line) for h = IO^Tk as 
a function of uj for U/ttA = 3.0 compared with the NRG results (dashed line). 

this is because a discrete spectrum had to be used, and the broadening, which is apphed on 
a logarithmic scale, is much larger at higher values of uj. 

Finally, we check the extreme limit, h/T^ = 10^ (In/i/Tx = 6.93), where it can be seen 
in figure ^2 that the renormalised parameters revert to their bare values. Again from the 
results given in figure ITHl it can be seen that the two sets of results correspond. In this regime 
the peak is at Up = 2h, and all many-body effects have disappeared. 

The similarity of the parallel susceptibility to the susceptibilities in the absence of a 
magnetic field suggests that the this response function should satisfy the Korringa-Shiba 
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FIG. 19: Results for ImXs,||(w, /i)/7rx^ ||(0, /i) calculated using the RPT equations for a range of 
values of the magnetic field h. 

relation. We test this out by plotting the imaginary value of h) divided by c<J7rx5_(0, h) 

in figure (Till for U/ti/S. = 3.0 and h/T^ = 0.1,0.4, 1 (the parallel suceptibility is very much 
suppressed for larger values of h). We see the the Korringa-Shiba relation does apply and 
with increasing magnetic field the width of the central peak broadens, reflecting the suppres- 
sion of the Kondo resonance with the increasing values of the magnetic field. The increase 
of area under the peak with increase magnetic field, from equation (|T^ . can be related to 
the corresponding decrease in the static parallel susceptibility. 



IV. CONCLUSIONS 



We have taken into account processes involving repeated quasiparticle scattering to cal- 
culate the dynamical charge and spin susceptibilities at zero temperature. As we use the 
exact irreducible vertices at = 0, we know from Fermi liquid theory that the results will 
be asymptotically exact as —>■ 0, and that the Korringa-Shiba relation will be satisfied. 
We have also shown that the total spectral weight is given correctly. However, it is not 
self-evident that taking into account these processes alone in the RPT will constitute a good 
approximation over all the relevant the energy scales. By comparing the results with those 
obtained by a direct evaluation of the response functions using the NRG we have established 
that the results of the RPT calculations give a remarkably accurate description of the dy- 
namic spin susceptibilities for the Anderson model with U > over all the relevant energy 



scales and arbitrary values of the external magnetic field. The earlier work of Kuramoto 
and Miyake !12] came to the same conclusion, though their derivation is restricted to the 
Kondo limit, and with no magnetic field. 

Whereas the spin excitations are enhanced and pushed to lower frequencies for large 
values of U, the low energy charge excitations are suppressed and pushed to higher energies. 
The RPT results in this case can only accurately describe the low energy suppression, as the 
high energy peaks in the charge spectrum are pushed to energies of the order uj ~ U/2, where 
the low energy renormalised parameters used are not appropriate. It is the spin fluctations, 
however, which are the important fluctuations in this regime, which are enhanced and pushed 
to lower energies. The energy scale of these fluctuations in the large U regime is the Kondo 
temperature Tk, which is the renormalised energy scale, and over this scale the RPT is an 
excellent approximation. For the negative U Anderson model the opposite will be the case. 
For large absolute values of U the spin fluctuations are suppressed and the charge fluctuations 
enhanced and driven to low energies. In this case it will be the charge fluctuations that are 
physically the more important ones, and for these the RPT will provide a good description. 

In calculating the quasiparticle contribution to the low energy spectral density of the 
single particle Green's function G'o-(co'), it has to be multiplied by the quasiparticle weight 
factor z (wave function renormalisation factor). At first sight it might seem surprising that 
the quasiparticle contribution to the dynamic susceptibility does not have to be multiplied 
by a factor z'^. The reason why such a factor is not required is because there is also a 
renormalisation of the vertex A{uJi : u) which couples the system to the external magnetic 
field. It can be shown via a Ward identity that this vertex evaluated at tui = = is equal 
to 1/z [2^ dJid, as a consequence, on differentiating twice with respect to H brings in a 
factor 1/2;^ which cancels factors of z from product of the full Green's functions. As a result 
the z does not appear in calculating the dynamical response functions from the RPT. 

The approach used here can be generalised to lattice models, such as the Hubbard model 
and Hubbard-Holstein models, within the dynamical mean field theory, and these calcula- 
tions are in progress. 
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V. APPENDIX 

The free local quasiparticle propagator Gcr{uj) for the symmetric Anderson model, with 
a magnetic field H acting at the impurity site, is given by 

GaA^) = \ , ~ (29) 

u + zAsgn(a;j + an 

in the zero temperature perturbation theory formalism There are two indepndent quasipar- 
ticle pair propagators defined by 

UZiu, h)^— G^iu + u', h)G,,{u\ h)dw', (30) 

ZTT J-oo 

and 

h)^— G,{u - J, K)G,,{J, K)dw'. (31) 

^ in J-oo 

These integrals can be evaluated analytically and the results for a; > are 

fllliuj, h) = - ^ - for a; = 0, 

' In — ^ + In — ^ U for u; 0. (32) 



and 



7ru;(a; + 2iA) { \ i/!i,-h J \ iA + h 

nl]iu,h) = .^-^ l^irif'A^] for uj = -2h, 

''^^ ' 7r{iA-h) 27rA \iA + hJ 



TT [uj + 2h + 2iA V iA + h J uj + 2h \ iA-h 

(33) 

The formulae for a; < are obtained by replacing i with —i in the above. For particle-hole 
symmetry, we have n^J(u;, h) = -tll\{uj, h), and nl\{uj, h) = nf^(u;, /i) = tll\{uj, -h) 
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